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Abstract. We use a discrete multiscale analysis to study the asymptotic integrability of 
differential-difference equations. In particular, we show that multiscale perturbation techniques 
provide an analytic tool to derive necessary integrability conditions for two well-known dis- 
cretizations of the nonlinear Schrodinger equation. 

1. Introduction 

The nonlinear Schrodinger (NLS) equation, 

id t f + d xx f = a\f\ 2 f, f = f(x,t), a = ±l, (1) 

is a universal nonlinear integrable partial differential equation for models with weak nonlinear 
effects. Here x is the spatial variable and t is the time, while d denotes differentiation with 
respect to its subscript. It has been central for almost fourty years in many different scientific 
areas and it appears in several physical contexts, see for instance [4,5,28]. 

In their pioneering work [32] Zakharov and Shabat proved its integrability by solving its 
associated spectral problem. From the integrability of Eq. ([T]) it follows the existence of infinitely 
many symmetries and conservation laws, and the solvability of its associated Cauchy problem. In 
correspondence with its symmetries one finds an infinite number of exact solutions, the solitons, 
which, up to a phase, emerge unperturbed from the interaction among themselves. 

On the other hand, also the problem of the discretization of the NLS equation has been the 
subject of an intensive research. In literature, one may find a few discretizations of the NLS 
equation. An integrable differential-difference equation discretizing Eq. JT]) has been found by 
Ablowitz and Ladik [2]. It reads 

• r» p , fn+1 ~ 2fn fn—1 i r \2 fn+1 + fn— 1 , n \ 

idtfn + = a \M 2 ' V) 

where n £ Z and h is a real parameter related to the space-discretization. Eq. ^fy admits a 
Lax pair and it has an infinite number of generalized symmetries and local conservation laws, 
which provide an infinite number of explicit soliton solutions [4]. As one easily see, in the limit 
as h — > 0, Eq. goes into the NLS equation (TjQ). 

l 
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From the point of view of the applications, the most relevant differential-difference NLS 
equation is given by 

• q £ i fn+l ~ 2f n + fn— 1 i r i2 i /o\ 

mfn + ^2 = fn ■ ( 3 ) 

By introducing the parameter s = 0, 1, the discrete NLS equations (|2][3]) may be combined in 
the equation: 

.o , (/„+! -2/ n + /„_i) (1 -S(T/l 2 |/ n | 2 ) 2 

idt/n + =o\fn\ fn- (4) 

The case s = 1 corresponds to Eq. ([2]) while the case s = gives Eq. ([3|). 

Eq. ([3]) is one of the most studied lattice models (see for instance [1,4,9,13-15] and references 
therein). Its study has a long and fascinating history, beginning in the 50's of the last century in 
solid state physics with the Holstein's model for polaron motion in molecular crystals [19]. Later 
on it appears in biophysics with the Davydov's model for energy transport in biomolecules [27]. 
Among the many recent applications of Eq. ([3]) let us just mention the theory of Bose-Einstein 
condensates in optical lattices [1] and semiconductors [15]. Its continuous limit goes again 
into the integrable NLS equation (pQ). The discrete NLS equation ([3]) possesses exact discrete 
breathers solutions [14] and just a few number of conserved quantities and symmetries are 
known. Numerical schemes have been used to exhibit its chaotic behavior [3]. A proof of its 
non-integrability, based on multiscale techniques, has been recently presented by the authors 
in [25]. 

Multiscale analysis [29, 30] is an important perturbative method for finding approximate 
solutions to many physical problems by reducing a given partial differential equation to a simpler 
equation, which can be integrable [8]. Multiscale expansions are structurally strong and can be 
applied to both integrable and non-integrable systems. Zakharov and Kuznetsov [31] have shown 
that, starting from an integrable partial differential equation and performing a proper multiscale 
expansion, one may obtain other integrable systems. In particular, they showed that the slow- 
varying amplitude of a dispersive wave solution of Eq. (fT]) satisfies the Korteweg-de Vries (KdV) 
equation, and viceversa. 

Let us give a sketch of their derivation, showing how to obtain a KdV equation as the 
lowest order of the multiscale expansion of the NLS equation (TIJ with a = 1, the so called 
repulsive NLS equation. One separates the complex field / in its amplitude and phase, f(x, t) = 
[v(x, t)/2] 1 / 2 exp[i(j)(x, £)], and rewrite the NLS equation as the system of real partial differential 
equations 

dm + d x (utp) = 0, (5) 

dttp + <pd xV > + d x u = ^dxv-^dy/ 2 , (6) 

where ip = d x (j). For long waves and small perturbations around the equilibrium solution we can 
define the following formal perturbation expansions: 



u(x,t) = l + XWV.f), (7) 
1=1 

oo 

<P(x,t) = -i + j^-V'V.O, (8) 



i=l 



where x' = e(x — t) and t' = e 3 i are suitable slow- variables, since e is a small perturbation 
parameter. By inserting expansions ([7][H]) into Eqs. (J2K])> a direct computation shows that the 
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lowest nontrivial order of the perturbation, that is e , provides an evolution equation for the 
field v^-' w.r.t. the slow-time t': 

d t ,vM + * w^w - = 0, 

A o 

that is a KdV equation. 

In [8] Calogero and Eckhaus have used the multiscale technique, at its lowest nontrivial 
order, as a tool to give necessary conditions for the integrability of large classes of partial 
differential equations both in 1+1 and 2+1 dimensions. In particular it has been shown that the 
non-integrability of the resulting multiscale reduction is a consequence of the non-integrability 
of the ancestor system. The derivation of the higher order terms of multiscale expansions 
has been carried out by Degasperis, Manakov and Santini in [11] and Kodama and Mikhailov 
in [21]. In [12] Degasperis and Procesi introduced the notion of asymptotic integrability of 
order n by requiring that the multiscale expansion be verified up to a fixed order n of the 
perturbation parameter. An integrable partial differential equation, as the NLS equation ([1]), 
has an asymptotic integrability of infinite order. 

Recently, some attempts to extend this approach to discrete equations have been proposed 
[6,16-18,22-24]. In [22-24] a multiscale technique for dispersive Z 2 -lattice equations has been 
developed. The main idea of the method used in [22-24] is based on dilation transformations 
of discrete shift operators. Up to our knowledge this has been carried out for the first time by 
Jordan [20] . Let us illustrate the basic procedure in the case of a function f n : 7L — > C depending 
only on one discrete index. Let T n be the shift operator, T n f n = f n +i, and A n = T n — 1 be 
the difference operator of order one. The difference of order j in a new discrete variable n' is 
expressed in terms of an infinite number of differences on the lattice of variable n by means of 
the following formula [20]: 




where uj is the ratio of the increment in the lattice of variable n with respect to that of variable n' 
and the coefficients and S 3 k are the Stirling numbers of the first and second kind respectively. 

The Jordan formula ([9]) implies that a rescaling of a lattice variable gives rise to nonlocal 
results. Therefore, to avoid the presence of infinite sums one needs to truncate the series ([9]), 
namely to introduce a slow-varying condition: 

A P „ +1 /n = 0, (10) 

p being a positive integer. A more general slow-varying condition has been recently introduced 
in [26]. 

Indeed, in [22-24], the multiscale analysis has been performed taking into account condition 
(fTUl) . As a consequence, the reduced discrete equations turned out to be non-integrable even if the 
ancestor equation was integrable. However, as shown in [16-18], if p = oo the reduced equations 
become formally continuous and their integrability may be properly preserved by the discrete 
multiscale procedure. In this way multiscale techniques easily fit with both difference-difference 
and differential-difference equations. The results contained in [16-18] confirm a discrete analog 
of the Zakharov-Kuznetsov claim [31]: "if a nonlinear dispersive discrete equation is integrable 
then its lowest order multiscale reduction is an integrable NLS equation" . 

In this paper we present the multiscale perturbation analysis of Eqs. d2H3|), thus extending 
to the discrete setting the approach used in [11,12,21,31]. The derivation of the higher order 
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terms in the perturbation expansion will enable us to provide an analytic evidence of the non- 
integrability of Eq. (J3j) - In fact, even if its lowest order reduction is an integrable KdV-type 
equation, the higher orders reductions exhibit non-integrable behaviors (see also our recent 
letter [25] where no details were presented). On the contrary, the same calculations for the case 
of Eq. ([2]) will show that the Ablowitz-Ladik discrete NLS equation satisfies all the integrability 
conditions up to the same order considered in the non-integrable case. This is an indication of 
its asymptotic integrability of finite order but not a proof of its integrability as for it we should 
go up to infinite order. 

The paper is organized as follows. Section [2] is devoted to the presentation of some technical 
details and basic formulas for the multiscale analysis of Eqs. ([2j{3]) . The main results of the 
perturbation analysis will be given in Section El In the concluding Section O we discuss further 
perspectives of this approach. 

2. Basic formulas for multiscale analysis of discrete NLS equations 

As for the continuous NLS equation, also for the discrete NLS equation fl3J) we introduce 
amplitude and phase of the function f n (t), namely f n (t) = [ z/ n(i)] 1 ^ 2 exp[i</> n (t)]. Therefore the 
discrete NLS equation may be written as the following nonlinear system of real differential- 
difference equations (s = 1 for ([2]) and s = for ([3])): 

d t v n = (sau n - [yjv n v n+l sin(^ n+ i - n ) + ^/zT^v^I sin (</>„_ i - n )] , (11) 



1 

^2 + (s - 2)au n 



COs((/> n +l - 4> n ) + J 1 ^—^- COs(^ n _l 



• (12) 



By analogy with the continuous case, see Eqs. (015]), the real fields v n (t) and <j) n (t) are 
expanded around the constant solution f n (t) = exp (— iat) in the following way: 



oo 



U n (t) = l + Y,£ 2l ^(K,{tm} m >l), (13) 
i=l 

oo 

</> n (t) = -rt + ^ £ 2,4 ^( K ,{U m > 1 ) 1 (14) 
i=l 

where e, with < e <C 1, is the perturbation parameter. The fields i/W and 0® in Eqs. (113H14I) 
depend on the slow-space variable k = e£n, ( £ R, and the slow-time variables t m = e 2m_1 t, 
m > 1. The free parameter £ will be fixed later so as to obtain a suitable continuous limit. 

In general, given a function u n {t) = v(n, {t m } m >i) we expand u n ±i(t) and dtu n (t) in terms of 
the slow variables k and {i m } m >i (see [16,17] for further details). Let T n be the shift operator 
defined by T^u n = u n ±\. Then we have: 

u n ±i = (T^Y v(k, {t m }m>i) = ^2 ~^~ V ( K ' {Wm>l) , (15) 

with 



i=0 



f-lV" 1 

^ = EH^ A - A* k = (T k -1)\ (16) 



and 



d t u n = J2 e 2i - l d ti v(K, {tmU^). (17) 



i=l 
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If u n is a slow- varying function of order p, see Eq. (jlOp , we can truncate the infinite series in 
Eq. (|16p . In such a case the 5 K -operators reduce to polynomials in the A K -operators of order at 
most p. Hereafter we shall assume p = oo and the <5 K -operators are formal differential operators. 

Taking into account the expansions (|13til4p and Eqs. (I15p and (I17p we have the following 
formulas for the shifts of the functions v n {t) and (f) n (t), 

oo UM . .. ,j_2i 
u n±1 = l+E e 'E (j _2i)! V®Mtm}m>l), (18) 



j=1 i=l 



and their time derivatives, 



~ CT * + E ei E + 1)! * w MW*gi), (19) 



oo j-1 

d t v n = J^E^^^'W^i), (20) 

i=2 i=l 

oo j 

a t n = -a + XI^'E^'-^C/e.^mW). (21) 

j=l i=l 

3. Main results 

The multiscale analysis of the system of real differential-difference equations (Illtil2p is carried 
out by inserting the formal expansions (|13til4p and (I18ti2ip into Eqs. (|llH12p and by requiring 
that the resulting equations be satined at all orders in e. 

The lowest non-trivial order corresponds to e 2 . It gives 

From now on all results will be presented only for the phase functions -i * — 1? since any 
i > 1, is obtained from the even perturbation orders and may be expressed in terms of the 
0UJ ' s with j < i and their derivatives. 
At order e 3 we get 

(^,-cW>=0, C=± ^-f>' /2 . 

As c has to be real our multiscale analysis is performed only for a = 1. Note that < h < 1. 
Moreover we choose ( = h, so that c = ±(1 — s/i 2 ) 1 / 2 ; note that c remains finite as h — > 0. 
Therefore the resulting equation at this order is satisfied by </>M = ^'(a?, {t m } m >2) with x = 
k — ct\. Here we are assuming that the solution is asymptotically bounded. 

At order e 5 , the no-secular term condition implies (d 2 — c 2 <5 2 ) (f)^ = 0, so that (f>^ = 
(j)( 2 \x,{t m } m > 2 ). At this same order, the evolution equation for (jP-' w.r.t. the slow-time t 2 
reads 

t2 ^ = K 2 , K 2 [</> (1) ] = ai dl<j>W + a 2 [d x ^y, (22) 

with 

c 3 
«i = — [3 - (3s + l)/i 2 ], a 2 = sh 2 - -. 

Eq. (|22|) is a potential KdV equation and K 2 [<p^] denotes exactly the second flow of the 
integrable hierarchy associated with the potential KdV equation. A necessary condition for 



D. LEVI, M. PETRERA AND C. SCIMITERNA 



the integrability of the system ()lltil2j) is that its multiscale reductions provide the integrablc 
evolution equations (j > 3): 



(23) 



where C is the recursive operator associated with the KdV hierarchy, 



C[f(x)]=d 2 J(x) 



dufiu) 



and the /3j-'s are real coefficients to be fixed later. 

According to a general procedure for the multiscale analysis of partial differential equations 
[10-12] we now assign a formal degree to the x-derivatives of the functions (jP\ 



deg ( d. 



+ 2j - 1, 



>0, 



and define V n as the vector space spanned by the products of all derivatives d^.<p^ with total 

degree n. We denote by Vn C V n the subspace spanned by those products of derivatives d x (j)^ 
with j < r. 

After caring for secularities, the order e 7 yields (p^ = (p^ 3 \x, {t m } m >2) and the following 
non-homogeneous evolution equation for the field 4>^ w.r.t. the slow-time t2, depending on 
and its derivatives: 



d t2 4> {2) - ai d 3 x ^ - 2a 2 a^ (1) ^ 2 ) = -d t J^ + a 3 (d 2 ^) 2 + a A (d : 



where 



"3 



"5 



/t 2 [16/t 2 s-5(l + 3s)] + 7 
64 ' 

h 2 [16h 2 s -3(3 + s)] -3 
48 ' 



a 'i 



«6 



ch 2 (l + 7s) 



(24) 



12 



c[h 4 (15s + 1) + 30h 2 (s - 1) - 15] 
1920 



Substituting Eq. (|23|) with j = 3 into Eq. (|24|) and fixing 0$ = —u% in order to remove 
residual secularities, Eq. (|24|) reduces to the following evolution equation for the field (f)^ w.r.t. 
the slow-time t%: 

d t2 ^-K' 2 [^]^ = f^, (25) 
where K'- [<fi^] ip is the Frechet derivative of the flow Kj [4>^] along the direction ip, 



ar 



r=0 



In Eq. (I25h the forcing term is a well-defined element of Vr , dim'Pg 1 ^ = 3, namely a 
linear combination of three independent differential monomials (see Appendix, Eq. (|36p ). with 
known coefficients which are polynomial functions of h. 

Now the request for integrability of (|lltil2p implies the existence of the following evolution 
equation for the field <p^ w.r.t. the slow-time t^: 



d t3 ^ - K> ^ = /<*> 



(26) 
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where G 5 dim'Pg 1 ' , = 6 (see Appendix, Eq. (|37|)). so that the following compatibility 



condition must hold: 



{d ta -K' 3 [<pW]}f^ = {d t2 -K!,[<pW]}f 



K*3) 



(27) 



Such a condition allows us to express the coefficients of the polynomial in terms of those 
of /(* 2 ) and it does not impose any further constraint on the coefficients of (see Appendix, 
Eq. (|38p ). As in our case this condition is satisfied, we conclude that the nonlinear system 
(Illtil2p has an asymptotic integr ability of order seven irrespective of the value of s. 

The next perturbation order, that is e 9 , gives rise to a bifurcation between the non-integrable 
(s = 0) and the integrable case (s = 1). For the sake of clarity we will study separately the two 
cases. 

• The case s = 0. At order e 9 , the resulting equations provide the evolution of the field 
0( 3 ) w.r.t. the slow-time t 2 . This is given by an integro-differential equation. To reduce 
it to a purely differential equation we introduce the fields <pj> = d x (f>^\ Taking care of 
secularities and taking into account that cfP-> evolves w.r.t. the slow-time t± according 
to Eq. (|23D with j = 4, we get <jW = 4>^{x, {t m } m > 2 ) and 



(i) 



(3) 



(28) 



where H'- [</?^] ip is the Frechet derivative along t/j of the j-th. KdV flow Hj [cp^ 



d x Kj [<p^~\- Here is a known element of the space Vg , dimT-g = 14 (see Ap- 
pendix, Eq. ([4*0]) ). The evolution equation of w.r.t. the slow-time t% takes the 
form 



9t 3 <4> 



(3) 



H' 3 



(i) 



(3) 



9 



(*s) 



(29) 



where the coefficients of g( t3 ' G V 
compatibility condition 



{?, dim^i? 



31, are determined by requiring the 



{dt 3 



H' 3 



(i) 



}/ 2) = { 



t-2 



H' 



(i) 



(30) 



Eq. (|3U|) is a necessary condition for the integr ability of the system (|1HI12|) with s = 0. 
In this case only nine out of the fourteen coefficients of g^ are independent. Thus 
we have five integrability conditions (see Appendix for further details). It turns out 
that the obtained constraints on the polynomial g^ are not satisfied by the coefficients 
computed in Eq. (|28p . Therefore the system (|!Hll2p with s = 0, namely the discrete 
NLS equation (f3j), does not fulfil the necessary conditions assuring its integrability. 

The case s = 1. In this case the resulting equations are purely differential and one can 
remain within the potential KdV hierarchy. Taking care of secularities and taking into 
account that (p^ evolves w.r.t. the slow-time t± according to Eq. (I23p with j = 4, we 

get = 4>( 4) (x,{t m } m > 2 ) and 



^(1)1 0(3) = h (t 2 ) 



where is a known element of the space Vf\ dim 'P. 
of cj)^ w.r.t. the slow-time t% takes the form 



(2) 



(31) 

11. The evolution equation 



d t3 <f>® - K ' 



(i) 



4® = , 



(32) 
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where the coefficients of G 'Pwi dim 7^ 
compatibility condition 



(2) 



24, are determined by requiring the 



d t3 - K!, 



d t2 - k', 



0«] } hM . 



(33) 



In such a case it turns out that all the constraints imposed by (j33[) on the eleven 
coefficients of the polynomial are satisfied by the coefficients computed in Eq. (|31|) 
(see Appendix for further details). This proves that the system (1 1 111 1 2 h with s = 1, 
namely the discrete NLS equation ([2j), has an asymptotic integrability of order nine. 
Actually, since the discrete NLS equation ([2]) is known to be integrable, its asymptotic 
integrability should be of order infinite. 

The above results may be summarized in the following Proposition. 
Proposition. The nonlinear differential- difference equation 

fn+1 — 2fn + fn-1 



id t fn + 



2h 2 



\fn\ fn 



is non-integrable. In particular, its multiscale reduction, carried out by using the formal expan- 
sions U3\\ 14\ ) and \18W21\) . shows that it has an asymptotic integrability of order seven. 
The differential- difference Ablowitz-Ladik equation 

fn+1 — 2fn + fn-1 _ , j. ,2 fn+1 + fn-1 



idtfn 



\fn\ 



2h 2 2 

has an asymptotic integrability of order nine. (Actually its asymptotic integrability should be of 
order infinite since it is known to be integrable.) 



4. Concluding remarks 

The present paper has been devoted to the derivation of higher order terms of the multiscale 
perturbation of discrete NLS equations around the constant equilibrium solution. This enabled 
us to study the asymptotic integrability of Eqs. (J2H3J) , thus proving that the discrete NLS 
equation (|3|) is non-integrable. Such a result has been already established in [25], but a detailed 
presentation of the integrability conditions appears for the first time in the present paper. 
Moreover, we have also investigated the asymptotic integrability of the Ablowitz-Ladik discrete 
NLS equation ©. 

We notice that the obtained results can be also used to construct approximate solutions of the 
discrete NLS equations (J2}|3]) . They will be expressed in terms of the solutions of the continuous 
equations belonging to the KdV and potential KdV hierarchies. More precisely, the solutions of 
the lowest order term of the multiscale expansion of (J2]{3j) will be expressed in terms of a soliton 
solution of the potential KdV equation. 

It is worthwhile to notice that the presented discrete multiscale technique fits with both 
differential-difference and difference-difference equations. Therefore it can be used to investigate 
the asymptotic integrability of a large class of discrete dynamical systems. The method turns out 
to be a useful analytic tool whenever one has to deal with a discrete equation whose integrability 
is not established yet. 

As a future work, we plan to investigate the asymptotic integrability of the following class of 
differential-difference equations [7]: 

o f fn+1 ~ 2f n + fn—1 f , f ft \ 

°tjn H TTTn — In + 9{In-l, In, In+l), 
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where g is a homogeneous polynomial of degree three. The importance of the above class of 
equations lies in the fact that it admits translationally invariant discrete kinks solutions. 



Appendix: The integrability conditions for the potential KdV and KdV 

hierarchies 

This Appendix is devoted to the presentation of the integrability conditions for the potential 
KdV and KdV hierarchies we used in our derivation. In particular, the following formulas have 
been used to obtain the results contained in Section [3l Thus we shall use the same notation. 

The potential KdV hierarchy. The integrable hierarchy of the potential KdV equation 
is given in Eq. (|23f) . The quantities Kj[cp^] and their corresponding linearizations K'j[(p^}ip, 
for j = 2 and j = 3, read (here d = d x ) 



K 2 



(i) 



aid 3 </> (1) + "2 



i) 



5a 2 



6a\ 



^ (V 1 )) 3 + (sy 1 )) 2 + 2d0( 1 W 1) 



(34) 
(35) 



and 



K> 2 



K' 3 



10a 2 
3ai 



dct>Wd 3 ip + d 2 <t) (1) d 2 i) + ^ (a 2 (1) ) dip + sP<f)Wd<ip 



where a%, a 2 , fh are real coefficients (in our case they are polynomial functions of the parameter 
h). 

The non-homogeneous terms 6 'Pq 1 ^, £ ^8 j gi ven i n Eqs. (|25l) and (|26|) respec- 
tively, are 



= a 1 (^ 1 )) 3 + a 2 ^ 1 )5 3 0( 1 )+a 3 (9V (1) ) 2 , 

/(*3) = (d 2 ^ 1 A\b 2 d<p^d 5 <p^ + b 3 d 2 ^dU {l) 



+ b 4 



i) 



+ 65 



(36) 



(37) 



The combatibility condition ([27]) implies the following algebraic relations between the coeffi- 
cients 01,02,03 and 6i,...,&6 : 

9af 61 = 5/3 3 [9aiai + 2 (a 2 + 3a 3 ) a 2 ] , 
3«i6 2 = 5/3 3 a 2 , 
3ai6 3 = 5/3 3 (a 2 + 2a 3 ) , 
54af 64 = 5/3 3 a 2 (27ai«i — a 2 a 2 ) , 
9af 65 = 5/3 3 (9aici!i + 5a 2 a 2 ) , 
k 3ai6 6 = 5/3 3 (a 2 + a 3 ) . 



(38) 



The system (|38|) allows us to express the b{S as functions of the a,'s without requiring any 
constraints on the latter ones. This means that the compatibility condition (|27[) is satisfied for 
any ai,a 2 ,a 3 provided that ([3"8]) is fulfilled. 
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The non-homogeneous terms G G "Piq > defined in Eqs. (fBTj) and (f3"2"j) re- 

(2) (2) 

spectively, are quite long, since dim'Pg ' = 11 and dimP]^ = 24. In order to present the 
integrability conditions imposed by the compatibility equation (f3"3"j) it is sufficient to write down 
only the expression of h^ t2 \ It reads 



ci 

+ C 6 



i)) 2 ay 1 ) + C6 ^(i)^ 4 + ^5^(1)^0(2) + C8 ^ ( i) ^ 



+ c 9 a 3 0W^ 2 )+c lo (d^f d<f>® 



en 



2) 



(39) 



The compatibility condition (|33|) allows one to express the twenty-four coefficients of M* 3 ) in 
terms of those of hS t2 > (these algebraic relations are easy to derive but rather cumbersome and 
we do not present them here) and the following three algebraic constraints involving the coeffi- 
cients ci, en and the coefficients oi, a%, 03, a±, a 2 previously defined. The obtained necessary 
integrability conditions read: 



C6 



[27oi (Q2 + 4a 3 ) ai - (37a| + 46a 2 a3 + 12a§) a 2 ] c n 
108af«2 

[(17o2 + 1803) a 2 — 27aiai] cs [(3a 2 — 803) a 2 — 3ai«x] c 9 

~T 77777 o ~t~ 



108a? 



36af 



(I8C2 - 24ci - 55c 3 ) a| (13c 5 - 3c 4 ) a 2 



54a? 



+ 



18ai 



C7 

cio 



a 2 
3aiCn 

a 2 



The KdV hierarchy. For the integrable hierarchy of the KdV equation we have 



and 



H 2 
H 3 



r 



(i) 



Ql dV 1 ) +2a 2 ¥> (1) 0¥> (1) , 



^2 



(i) 



(i) 



10a 2 
3ai 



ai V / 



dV 2) + 2a 2 



(i) 



p = a.\crip 



i) 



10a 2 
3ai 
2a 2 



Wd 3 P + 2d<pV>& P + ( 2aV (1) + - (V 1} H 5p 



+ , ^(1)^(1)+ S3 (1) 



The above expressions are obtained by differentiating with respect to x the corresponding 
expressions given in Eqs. (13411350 and setting = d x (f)^ l \ p = d x i(). 

The non-homogeneous terms G Pg\ 9^ £ ^11 > defined in Eqs. (|28|) and (|29|) re- 

(2) (2) 

spectively, are quite long, since dim'Pg = 14 and dimP]^ = 31. In order to present the 
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integrability conditions imposed by the compatibility equation (|30p it is sufficient to write down 
only the expression of g^" 2 ' . It reads 

g w = ^ay^V 1 ) + dzdvW&tpW + d 3 <p^d 5 v w + d4 (V 1 )) 3 

+ d 9 d<pWd 2 <pM + d 10 &<p<- 1 W 2) + W 2 W 1} + d 12 (tpwydtp® 

+ (WV 2) <M 1) + d 14 v (2) d<p (2) . (40) 

The compatibility condition (j30j) allows to express the thirty-one coefficients of terms 
of those of g( t2 ' (these algebraic relations are easy to derive but rather cumbersome and we do 
not present them here) and the following five integrability conditions involving the coefficients 
d%, du and the coefficients a±, a 2 , «3, ot\, a 2 previously defined: 

[9ai (12a 3 + 5a 2 ) a\ - (45a 2 + 88a 2 a 3 + 12a|) a 2 ] du 



d 7 



5Aa\a 2 

[(3a 2 — 803) a 2 — 3ai«i] dxo 2 [(21a3 + 4a 2 ) a 2 — 9aia±] dg 

+ 9^2 + ?f^x 

+ {9d 5 + 8d 6 - 24d 4 ) a 2 2 (12d\ - 30d 2 + 85d 3 ) a\ 



9«i 27aj 



dx 



a 2 d 
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2a 2 

a 2 du 



du = d w -d 9 + 
di 2 ■■ 



2a 2 

3o.i cZi4 



d 



13 



2a 2 ' 
3oi<ii4 
a 2 
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